Tailoring the frictional properties of granular media 
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A method of modifying the roughness of soda- lime glass spheres is presented, with the purpose 
of tuning inter-particle friction. The effect of chemical etching on the surface topography and the 
bulk frictional properties of grains is systematically investigated. The surface roughness of the 
grains is measured using white light interferometry and characterised by the lateral and vertical 
roughness length scales. The underwater angle of repose is measured to characterise the bulk 
frictional behaviour. We observe that the coefficient of friction depends on the vertical roughness 
length scale. 

PACS numbers: 45.70.-n 81.65.-b 



I. INTRODUCTION 

Experiments have confirmed that friction plays an im- 
portant role in the physics of static granular materials. 
It was demonstrated that the random loose packing vol- 
ume fraction decreases with increasing inter-particle fric- 
tion coefficient 1,2 . The role of friction on force chains 
was investigated in reference^. If static granular media is 
described by a statistical mechanics approach^, a "con- 
figurational temperature" can be defined, which depends 
on friction 5 . 

Friction plays a vital role in granular dynamics, too; 
for instance in the case of granular material undergo- 
ing shear——. The effect of friction in granular mate- 
rials under shear in a linear, split-bottomed shear cell 
filled with layers of high-friction and low-friction grains 
has been investigated^. Here, it was shown that the 
position of the shear band depends on the layering of 
the two types of grains. It has been shown experimen- 
tally that friction can also play an important role in 
granular segregation^^, for example in the sharp tran- 
sition from the reverse brazil-nut effect to the brazil- 
nut effect 10 , but other work shows that this is not al- 
ways the caseii. The dynamical properties of avalanches 
were shown to depend on inter-particle friction, and two 
friction-dependent scenarios for granular avalanches on 
an incline were identified— - —. 

However, none of the above work was made with par- 
ticles of the same material and shape but more than two 
coefficients of friction. In this study, it is shown that it is 
possible to tune the topographic properties of soda-lime 
glass grains, making inter-particle friction a control pa- 
rameter for some types of beads and thus allowing the 
systematic quantitative study of the principles underly- 
ing granular phenomena such as those mentioned above. 
Figure [1] shows the surfaces of grains when they are un- 
treated and after etching with the two protocols we de- 
scribe below. 

The article is structured as follows: in section [III the 
chemical etching protocols are described. The analysis 
of the grain roughness follows in section IIII[ In section 




FIG. 1: SEM images of the surfaces of MoSci glass spheres 
of diameter 146 /mi ± 19 /xm (small Mo-Sci). Left: smoothed 
by NaOH. Centre: untreated. Right: etched for four minutes. 
The white bars are 10 fim long. 



HVl we present the bulk frictional properties of the grains 
measured by the underwater angle of repose and in ap- 
pendix [C] a bulk measure of roughness is presented. 



II. ETCHING GLASS GRAINS 

Usually, hydrofluoric acid (HF) is used to etch glass^i^. 
However, HF alone is unpleasant to work with and is 
harmful to the environment^. Furthermore, our experi- 
ments with dilute HF showed highly inhomogeneous etch- 
ing unconducive to our purpose, as demonstrated in fig- 
ure [2j Ammonium bifluoride showed better etching prop- 
erties (figure [1]) as well as better handling: an aqueous 
solution of ammonium bifluoride will contain only small 
amounts of free HF, according to the set of equilibrium 
equations [TJ 
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FIG. 2: Optical micrograph of a glass sphere etched in HF. 
The white bar is 50 /im long. Instead of homogeneous surface 
roughening, large craters are produced. 



NH 4 FHF ^ NH+ + FHF" ^ NH+ + HF + F" (1) 
3HF + 2H 2 ^ FHF" + F" + 2H 3 + 

A protocol was developed using ammonium bifiuoride 
which optimises the flow properties of the etch solution 
whilst releasing little hydrofluoric acid and etching as 
homogeneously as possible 17 . 

A. Roughening protocol 

This is a hydrofluoric acid etching procedure, so the 
utmost care must be observed during all steps in order 
to avoid chemical burns. Since the procedure is one to 
etch glass, plastic lab ware should used. 

1. 14 parts by mass of liquid glycerine, C3H8O3, and 
six by mass of distilled water are heated in a water 
bath to 90°C. 

2. Six parts by mass of solid ammonium hydrogen- 
difluoride, NH4HF2, and one part by mass of solid 
granulate iron (III) chloride- hexahydr ate, FeCl3 • 
6H2O (a catalyst), are added. 

3. An over-saturated solution is produced by mixing 
vigorously at 90° C. 

4. The solution is cooled under agitation until it is 
cool enough to be placed in the centrifuge. 

5. The over-saturated solution is centrifuged until the 
clear solution is separated from the undissolved so- 
lute (about 15 minutes). 

6. The solution is added to the glass grains and shaken 
vigorously. After the desired etching time has 
passed, the solution is poured out and the grains 
thoroughly washed with distilled water. 



A separate etching procedure to make the grains 
smoother was also employed. This method was published 
by Schellenberger and Logan^, To smoothen the grains, 
they are immersed in 12.5 M sodium hydroxide, NaOH, 
for 30 minutes, under occasional agitation. 

We use soda-lime glass grains from two different man- 
ufacturers and in two different sizes. These are presented 
in the table below. 



Name 


Diameter (/xm) 


STD (/xm) 


Symbol 


Small Mo-Sci^ 


146 


19 


A 


Cataphote^ 


221 


40 


□ 


Large Mo-Sci 


240 


20 


• 



TABLE I: (Colour online) The types of grains used in this 
study. The particle diameter and standard deviation (STD) 
was measured using a Camsizer, as outlined below. 

We perform seven treatments on each type of glass 
grain: smoothed for 30 minutes in NaOH, as supplied 
by the manufacturer ( "unetched" ) and roughened for 30 
seconds and one, two, three and four minutes. 

The effect of etching on particle size was measured us- 
ing a Retsch Technology Camsizer. The particle size dis- 
tribution of 200,000 grains of each batch was measured. 
The reduction in grain diameter after four minutes' etch- 
ing with the roughening protocol is very small: between 
unresolvable within experimental error for the large Mo- 
Sci to a reduction of 3.5 fim ± 1.5/im (a reduction of 
2.4%) for the small Mo-Sci grains. The smoothing pro- 
tocol does not bring about a resolvable reduction in di- 
ameter. The width of the size distributions also remains 
constant within experimental error. 



III. GRAIN ROUGHNESS 

The grain surface topography was measured us- 
ing white-light vertical scanning interferometry (WLI), 
(ZeMapper, Zemetrics Inc., Tucson, Arizona)^. WLI is 
an optical surface scanning method that is able to pro- 
vide nm height resolution over a large field of view. This 
enables the characterisation and quantification of surface 
topography variations of, for example, natural grains, re- 
acted crystal surfaces, and mineral aggreg ates^. The 
principle of operation is shown in figure [3j After pass- 
ing through an external beam splitter, the white light 
passes through a temperature-stabilised Mirau objective 
(magnification: 160 x). The objective is equipped with 
a second, internal, beam splitter. One beam is reflected 
by a stationary reference mirror, the other by the sam- 
ple surface. The two light beams differ in optical path 
length as a function of the distance to the sample sur- 
face. The resulting interference pattern is captured by a 
high-resolution charge-coupled device (CCD). The focal 
plane of the interference pattern is scanned through the 
sample. From the collected interference pattern data, a 
height map is calculated. The vertical resolution of the 
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data set is <1 nm; the maximum field of view applied in 
this study is 80 /im x 80 /im (figures [3] and 3]) . 

For the WLI measurements, glass grains are bonded to 
glass microscope slides with ultra- violet-curable adhesive 
(Norland Optical Adhesive 61) and sputter-coated with 
40 nm gold. 



CCD 



white light source / \ 




beam splitter 



Mirau objective 



piezo-electric 
sample holder 



FIG. 3: (Colour online) Principle of operation of a vertical 
scanning interferometric microscope. 





FIG. 4: Height map of the cap of a large Mo-Sci grain etched 
for four minutes. 



A. Roughness analysis 

Characterising roughness alone is a complex endeav- 
our, and there is a variety of tools available to do so^. 
We use the height difference correlation functional. How- 
ever, when measuring the roughness of spheres, one is 
presented with additional, unique problems. To obtain 
reliable statistics, a large field of view is necessary. How- 
ever, the curvature of the grain means that the measure- 
ment method must have a large enough dynamic range 



FIG. 5: (Colour online) The method employed to mea- 
sure roughness takes account of the spherical topology of the 
grain. The height difference correlation, p(L), between pairs 
of points r(6,(f)) and r' (0\ <//), separated by their great-circle 
distance, L, is calculated. R is the radius of the grain. 



to be able to measure nano-scale features over a height 
range of several /im. WLI proved to be better than AFM 
for this purpose. 

Once a height map, z(x,y), has been obtained, this 
may not be treated as a plane, since it is a projection of 
a sphere surface onto a plane, which can never be simul- 
taneously conformal and distance conserving^. It is also 
necessary to make an angle-of-sight correction to account 
for having measured height in the z-co-ordinate instead 
of in the radial direction. This correction is explained 
in appendix [Al Additionally, a batch of soda-lime glass 
grains will not consist only of perfectly spherical grains. 

For each batch of grains, ten to fifteen grains are anal- 
ysed. From the surface relief measurement, the typical 
surface feature size characterised by a vertical, £ ver t, and 
a lateral length scale, £/ a t, is measured. 

For a one-dimensional interface, £ ver t, 6at an d the 
Hurst exponent, H, are obtained from the (second-order) 
height difference correlation function^ 

C(L) = ((h(x + L)-h(x)) 2 ) 1/2 , (2) 

where h is the height of the interface at position x or 
x + L. For the surface of a rough sphere described in 



spherical polar co-ordinates, the height difference corre- 
lation becomes 

p(L) = ((r'(8',4>')-r(0,4>)) 2 ) 1/2 , (3) 

where the points r(0, and r'(0', (j) r ) are separated by 
the correlation distance L. On a sphere, L is the great- 
circle distance between the two points calculated using 
the spherical cosine law: 

L = R arccos(cos 6 cos 6' + sin 6 sin 0' cos(0' — </>)), (4) 

where R is the sphere radius. The angle brackets indi- 
cate an ensemble average for one L. 

The height difference correlation, p(L), on a double 
logarithmic scale, yields the vertical and lateral satura- 
tion length scales, ^ ve rt and £/ a t, for a given etch time. 
This is the position at which a cross-over between a 
power-law type behaviour and a saturation plateau oc- 
curs, and is illustrated in figure [6] 
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FIG. 6: (Colour online) The roughness, p(L), on a double 
logarithmic scale, yields the vertical and lateral saturation 
length scales, £ ve rt and £j a t, here for a Cataphote grain etched 
for 30 seconds. The slope of the fit through the first five data 
points is the Hurst exponent, H. 

The roughness length scales are obtained from the WLI 
micrographs of the spheres by the following steps: first, 
the WLI measurement is performed over a field-of-view 
of 80 pm x 80 pm (60 pm x 60 pm for the small Mo-Sci). 
The radius of the grain, R, is measured by optical microg- 
raphy. The data is cropped to 60 pm x 60 pm (30 pm 
x 30 pm) to remove the bad pixels at the edges, then 
processed to remove spurious data adjacent to known 
bad pixel sites. Bad pixels make up less than 1% of the 



cropped data and they are mostly at the edges where the 
effect of grain curvature is highest. Next, the measured 
R and an estimate of the x and y positions of the "north 
pole" are taken as start values for a hemispherical fit to 
the data. Using the position of the centre of the grain, 
xo,yo and zo, obtained from the hemispherical fit, the 
data is converted to spherical polar co-ordinates r, and 

<t>. 

Following this, p(L) is calculated, taking account of 
the angle-of-sight correction (appendix [A]). Three million 
randomly-chosen pairs of points are taken to calculate 
p(L) for each grain. The vertical length scale, ^ ver t^ is 
obtained by laying spline through p(L); the length L^ink 
at which the slope of this spline becomes zero is detected. 
p(L) is averaged for L > Lkinki giving <f ve rt- A linear fit 
through the first five data points gives the Hurst expo- 
nent, H. The cross-over of the linear fit and the plateau 

invert) gives & a £. 

Strong asphericity in the grain will lead to distortion 
of the roughness analysis, since our analysis assumes a 
spherical topology. To establish the quality of the data, 
radial averaging is performed (see appendix [Bj and data 
from aspherical grains is discarded. Of the fifteen grains 
measured per batch, on average two were rejected for 
asphericity. 



B. Roughness results 

The effect of etching on the topography of glass grains 
is shown in figures [71 and [5J where the mean vertical and 
lateral roughness length scales, £ ver t and £i at respectively, 
are plotted as a function of NH4FHF etch time. 

From figures [3 and El one sees that there are two types 
of roughening etching behaviour: a monotonic increase in 
feature width and depth with etch time for the small Mo- 
Sci grains, where £ ver t increases at an approximate rate 
of 0.17 pm per minute and £i at at a rate of 1.4 pm per 
minute; and a broadening behaviour for the two larger 
types of grains, whereby surface features become wider 
but not deeper as a function of etch time. It is not clear 
from our data if the change in etch mechanism is con- 
nected to the different radii or changes in the glass com- 
position. The smoothed grains show no appreciable dif- 
ference to the unetched grains, and are omitted from the 
plots. 

The Hurst exponent may also influence grain 
interactions^. However, as can be seen from figure [6j 
the power-law regime of p(L) exists over less than an 
order of magnitude, so caution must be exercised in in- 
terpreting the behaviour of the Hurst exponent, H, as a 
function of etch time (figure [9]) . We see < H < 1 for 
all 21 grain types. H shows a slight downwards trend 
with increasing etch time which is the same for all three 
grain sets, suggesting that the Hurst exponent can also 
be tuned by the etching protocols presented here. 

It is desirable to be able to measure surface roughness 
in bulk samples. A possible method for doing this, using 
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FIG. 7: (Colour online) The mean vertical roughness length plQ g . (Colour online) The mean lateraJ roughness length 

scale, Urt, as a function of NH 4 FHF etch time. In this and scal6) ^ ag a function of N h 4 FHF etch time. 

all subsequent plots, the error bars are the standard deviation; 

for each grain type ten to fifteen grains were measured. The 

linear fit through small Mo-Sci points and the means of t; ve rt 

are in the range of 30 s to 4 min. 



a carbonated soft drink, is demonstrated in appendix IU1 



IV. BULK FRICTIONAL PROPERTIES 

In this section, the angle of repose, a, is measured to 
characterise the bulk frictional behaviour; to eliminate 
electrostatic and capillary attraction between grains, the 
measurements are made under water. 



A. Angle of repose measurements 

The angle of repose of each sample is measured using 
the set-up described in figure [TUJ A cell 25 cm wide, 
30 cm high and 40 grain diameters deep is filled with 
water 28 . The water temperature is kept at 21 ±0.5°C for 
all experiments to keep the viscosity constant. An upper 
reservoir is filled with the granular sample. Then grains 
are allowed to flow to the bottom part of the cell at a rate 
of approximately 140 grains/s. After an avalanche has 
occurred, the flow of grains to the pile is stopped. The 
pile is photographed with a black-and-white CCD with 1 
Megapixel and 8-bit depth. Each photograph is binarised 
using a threshold that is chosen automatically to fall be- 
tween the peaks of its bimodal grey-level histogram. The 
top and bottom edges of the pile are detected and the 
location of the linear part of the pile is detected by lay- 
ing a spline through the line describing the top of the 
pile. Within this range, linear fits are made to the top 
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FIG. 9: (Colour online) The mean Hurst exponent, H, as a 
function of NH4FHF etch time. 



and bottom of the pile; the angle between them is the 
angle of repose, a. The measurement was repeated ten 
to twenty times per grain type. The resolution in a is 
±0.15° 
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FIG. 10: (Colour online) Angle of repose measurement 
method, (a) A sketch of the set-up: a cell [1] 25 cm wide, 
30 cm high and 40 grain diameters deep is filled with water. 
The reservoir [2] is filled with the granular sample. The fun- 
nel rod [3] is raised by a screw thread to open the hole enough 
to let grains flow to the bottom part of the cell [4]. After an 
avalanche has occurred, the rod is lowered enough to induce 
granular arching, thus stopping the flow of grains, (b) CCD 
image of the pile. Linear fits to the bottom and top edges of 
the pile give the angle of repose (see text). 



Results 
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FIG. 11: (Colour online) Underwater angle of repose as a 
function of roughening etch time. The smoothing protocol, 
NaOH, is plotted on the negative x-axis. The error bars are 
the standard deviation. The coefficient of friction is computed 
by \i — tana, where a is the angle of repose^. 



Figure [TT] shows a as a function of roughening etch 
time. For the small Mo-Sci grains, the angle of repose 
increases monotonically from about 24° to over 27° as 
a function of etch time. As in the roughness measure- 
ments, the smoothed grains do not show an appreciable 
difference to the unetched grains. 

The Cataphote and large Mo-Sci grains again show a 
different behaviour to the small grains. Again, their be- 
haviour is a step function: unetched and smoothed grains 
have an angle of repose about one degree less than grains 
etched in NH4FHF for 30 s or longer. The appearance of 
these trends is similar to £ ver t as a function of etch time 
shown in figure [71 However, the Cataphote grains, whose 
^ ve rt is always lower than the large Mo-Sci grains, have 
a consistently larger angle of repose. 

One reason for this discrepancy could be the different 
fractions of non-spherical grains present in the Cataphote 
and Mo-Sci samples. The average length-to- width ra- 
tio, l/b ave , of the unetched Cataphote and large Mo-Sci 
grains was measured using a Retsch Technology Cam- 
sizer (sample size of 7,000 grains). With l/b ave = 1.17, 
the Cataphote grains show a slightly greater average as- 
phericity than the large Mo-Sci grains, l/b ave = 1.12. It 
has been demonstrated that higher asphericity leads to 
a larger angle of repose^. Another reason for the dis- 
crepancy could be that the two larger types of grain also 
produce different packing densities, whish would be 
a natural consequence of their aspherity^. Though we 
did not measure <p in this set-up, we show that (j) has 
a strong influence on the angle of repose (appendix ID]). 
Again, for the larger grains, the smoothing protocol has 



no apparent influence on a. 

A valuable result of this work is to show how the 
frictional properties of grains depend on their (tailored) 
roughness. This is illustrated in figure H2l the angle of re- 
pose has a dependency on £ ver t, albeit a weak/noisy one. 
The angle of repose showed no appreciable correlation 
with £i at , again suggesting that it is the depth of sur- 
face features, not their width, that controls the friction 
between grains. However, there is not a simple, linear 
dependency. 
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FIG. 12: (Colour online) Angle of repose, a, as a function of 
correlation length scales (a) £ ver t and (b) t;i a t- 
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CONCLUSION 



In this study, a method of modifying the topographical 
properties of soda-lime glass grains was presented, with 
the purpose of tuning inter-particle friction. We pre- 
sented a roughening and a smoothing chemical etching 
procedure. The surface roughness of the grains was mea- 
sured using white light interferometry and characterised 
using the height difference correlation function. Two 
types of roughening etching behaviour were observed: a 
monotonic increase in feature width and depth with etch 
time for the small grains; and a broadening behaviour 
for the two larger types of grains, whereby surface fea- 
tures became wider but not deeper as a function of etch 
time. We concluded that, when changing the roughness 
of glass grains, the type of grain matters. With the three 
types of grains we used, we are unable to distinguish 
whether this is due to size or material composition. The 
smoothing protocol made little difference to the surface 
roughness. The underwater angle of repose, a, was mea- 
sured to characterise the bulk frictional behaviour. For 
the small grains, a increased monotonically as a function 
of etch time. The two larger types of grains showed a dif- 
ferent behaviour: unetched and smoothed grains had an 
angle of repose about one degree less than etched grains. 
In the appendix, we presented a suggestion for measuring 
grain roughness in bulk samples. 

The ammonium bifluoride etching technique presented 
here allows the experimentalist to modify the roughness 
of glass spheres. Grains with tunable roughness are the 
ideal work-horse for experimental studies of wet (or "hu- 
mid") granular systems^ or of hydrodynamic interac- 
tions between grains in viscous fluid a 34 i 35 . It is shown 
that changing the surface roughness of grains directly 
changes the way they behave in the bulk. Besides the 
systems already shown experimentally to depend on fric- 
tion, which were reviewed in the introduction, this has 
other important consequences for granular community in 
that it paves the way for comparison with simulations, 
where friction between grains is a typical control param- 
eter. Additionally, the ability to tailor glass grains will 
allow systematic study of the role of friction in the static 
and dynamic behaviour of granular material, for instance 
granular segregation and pattern formation, avalanching, 
and (shear)flow, packing and force chains. 



Appendix A: Angle-of-sight corrections 

When analysing the roughness of the glass spheres, we 
make an angle-of-sight correction to account for our hav- 
ing measured height in the z-co-ordinate instead of in 
the radial direction. This is shown in figure [131 We cor- 
rect the vertical height, Az(x, ?/), to a radial feature size, 
Ar(0,0), by 



Az(x,y) 




FIG. 13: (Colour online) When analysing the roughness of 
the glass spheres, we need to make an angle-of-sight correction 
to account for our having measured height in the z-co-ordinate 
instead of in the radial direction. R is the sphere radius and 
Az(x,y) = h(6, <j)) is the height of the grain profile measured 
vertically from above. The correct radial height of the feature 
is given by Ar — 

° J cos 



The relevance of this correction can be estimated for 
the worst-case scenario: a feature at the outermost corner 
of the region of interest on a small Mo-Sci grain. In this 
case, = 37°; thus omitting this correction results in an 
error of 26%. 



Appendix B: Radial averaging to test for aspherical 
grains 

Strong asphericity in the grain will lead to distortion 
of the roughness analysis, since our analysis assumes a 
spherical topology. In our case, this is particularly true of 
the two smoothest sets of grains. To check if the grain is 
spherical enough that the roughness data is not falsified, 
radial averaging of the relief measurement of the grain is 
performed: 



(Ar(0))U. 



(Bl) 



Ar 



cos# 



(Al) 



For a perfect sphere, both the radial mean and the 
standard deviation are zero for all 0. For rough grains, 
small deviations from zero are expected, for aspherical 
grains large deviations. Since our analysis depends on 
our ability to locate the kink in p(L), it is vital that the 
mean and the standard deviation are smaller than the 
length scales we wish to measure, ^ ver t and £i at . We set 
the criterion that, if the radial mean OR standard de- 
viation is greater than £ vert for any angle corresponding 
to a radial distance S < 5^ a t, we reject the grain from 
our analysis, on the grounds that its asphericity has an 
influence on its p(L) at the relevant length scales. This 
is illustrated in figures HH and [151 
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FIG. 14: (Colour online) Radial average and standard devia- 
tion of a large Mo-Sci grain etched for one minute. This grain 
is accepted: both lines remain within the range [—£vert—£vert] 
for L < b^iat- 



studied systematically by T. S. Coffey 3 -^. She concluded 
that the surface roughness of the sweets is one of the main 
causes of the reaction; namely, surface features provide 
nucleation sites for bubbles of carbon dioxide. We use 
this principle in a prototype experiment for measuring 
surface roughness in bulk granular samples using a car- 
bonated soft drink. 

The experimental method is as follows: 4 ml ± 0.05 ml 
of a freshly-opened carbonated soft drink (Coca-Cola) is 
poured into a measuring cylinder. Pouring is smooth to 
minimise premature nucleation of bubbles. 500 mg ±0.5 
mg of granular sample (small Mo-Sci) is added through 
a funnel. The maximum volume of cola, beads and foam 
is measured and the volume of the grains and the initial 
volume of cola subtracted. 

Figure [16] shows that there is a trend: the cola produces 
more foam when it comes into contact with the more 
strongly etched particles. Of course, not just surface 
roughness but also other effects, such as surface activ- 
ity, buoyancy and occupation density of nucleation sites 
play a vital role in the reaction 3 -^. It is interesting to 
speculate on how this method could be calibrated. 
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FIG. 15: (Colour online) Radial average and standard devi- 
ation of a different large Mo-Sci grain etched for one minute. 
This grain is rejected: the standard deviation exceeds the 
range [—£ V ert—€vert] for L < 5£j at . 
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FIG. 16: Maximum volume of cola foam as a function of etch 
time (small Mo-Sci). 



Of the fifteen grains measured per batch, on average 
two were rejected for asphericity. 



Appendix C: Using a carbonated soft drink to 
measure the roughness of grains 

The popular experiment in which a geyser is produced 
by dropping Mentos into a bottle of Diet Coke 36 was 



Appendix D: Angle of maximum stability as a 
function of packing fraction 

The angle of repose (pile slope after an avalanche, a) 
is dependent not only on the friction between grains but 
on their packing density. The same is true of the angle of 
maximum stability (pile slope just before an avalanche, 
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otmax)-, which proved to be the more accessible quantity 
in investigating this dependence. 

The density of the grain material is determined to 
within 0.45 /ig/cm 3 with a Micromeritics gas pyknome- 
ter. The grains are added to a cell of square cross-section 
4 cm x 4 cm which is filled with water. A packing of 
grains is (re)produced on a Ling Dynamic Systems V409 
shaker. The packing fraction is determined to within 
an error of ±9 x 10 -3 by measuring the height of the 
packing in the container. The container is transferred 
to a board which is hinged at one end and raised slowly 
at the other, under computer control, by a stepper mo- 
tor. A black-and-white CCD with 1 Megapixel and 8-bit 
depth is used to record the surface of the granular sample 
while the angle of the container is increased. The camera 
and an LED light source are mounted such that they do 
not move relative to the container, thus avoiding spurious 
results from shadows; the camera is triggered as the step- 
per motor starts. Image processing is used to determine 
&max with an accuracy of ±0.5° as follows: First, subse- 
quent images of the surface of the sample are subtracted. 
The subtracted images are binarised; the threshold for 
binarisation is taken as the mean grey value for a non- 
moving bed plus four standard deviations. Then, con- 
nected moving pixels are detected. A continuous moving 
area of more than one grain (more than 50 pixels with 
our optics and set-up) is our definition of an avalanche. 
The angle at which this occurs is the angle of maximum 
stability, a max . 



The smoothest and the roughest samples of the small 
Mo-Sci grains were tested. In both cases, a clear de- 
pendency of a max on the packing fraction, 0, was ob- 
served. The effect was strongest for the roughest grains: 
here, the angle increased from 29° to 41° in the range 
0.55 < 6 < 0.58. 



A systematic study of the underwater flow of frictional 
grains as a function of initial packing fraction was made 
by Pailha et al 38 . For fixed inclination angles, they ob- 
served that avalanching occurs more readily for lower 
packing fractions, a result with which our figure [T7] is 
consistent. 
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FIG. 17: (Colour online) angle of maximum stability, a m ax, 
as a function of packing fraction for two types of grains: small 
Mo-Sci etched for four minutes (purple closed symbols) and 
smoothed in NaOH (green open symbols). 
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